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ABSTRACT 

In order to develop an infection, diarrhogenic 
Escherichia coli has to pass through the stomach, 
where the pH can be as low as 1. Mechanisms that 
enable E. coli to survive in low pH are thus poten- 
tially relevant for pathogenicity. Four acid response 
systems involved in reducing the concentration of 
intracellular protons have been identified so far. 
However, it is still unclear to what extent the regu- 
lation of other important cellular functions may be 
required for survival in acid conditions. Here, we 
have combined molecular and phenotypic analysis 
of wild-type and mutant strains with computational 
network inference to identify molecular pathways 
underlying E. coli response to mild and strong acid 
conditions. The interpretative model we have 
developed led to the hypothesis that a complex 
transcriptional programme, dependent on the two- 
component system regulator OmpR and involving 
a switch between aerobic and anaerobic metabol- 
ism, may be key for survival. Experimental validation 
has shown that the OmpR is responsible for 
controlling a sizeable component of the transcrip- 
tional programme to acid exposure. Moreover, we 
found that a \ompR strain was unable to mount 
any transcriptional response to acid exposure and 
had one of the strongest acid sensitive phenotype 
observed. 

INTRODUCTION 

The acidic barrier of the stomach represents a strong chal- 
lenge for many pathogenic enterobacteria (1). Therefore, 
although the correlation between acid resistance and 
pathogenicity is not absolute (2), the ability of some 
Escherichia coli strains to survive exposure to strong 



acid conditions is potentially relevant for pathogenicity 
(3-5). For this reason, the molecular and physiological 
response to acid stress has been the subject of intense in- 
vestigation (5). Four acid stress response systems (ARs) 
that can protect E. coli from low pH are known to date 
(5-7). Three of these depend on the external supply of 
amino acids (glutamate, arginine and lysine) and have 
been proposed to share the same basic mechanism: reduc- 
tive decarboxylation of the amino acid (consuming a 
proton) followed by extrusion of the product from the 
cytoplasm by a dedicated antiporter that also imports 
the original amino acid (5). One of the acid response 
systems, AR1, is active in the absence of amino acids (8) 
and is based on the FoFl ATPase (6,9). At pH 2.5-3 the 
pump extrudes protons from the cytoplasm, with the con- 
sumption of ATP (7). The details of this system are poorly 
understood, but it is known to be regulated by both the 
sigma factor RpoS and the catabolite repressor protein, 
CRP (10). Also the other acid response systems are 
regulated by RpoS (10,11), which makes the acid stress 
response a growth phase dependent process (5,12,13). 
Other genes found in a 15-kb region around gadA 
termed the Acid Fitness Island (AFI) have also been 
shown to be required for acid resistance in E. coli 
(11,14-16). Among the regulators acting on the AR2 
genes are the regulator proteins YdeO, GadE, GadX, 
GadW and the two-component system EvgA/EvgS 
(5,17-20). In addition, the system is influenced by cross- 
talk with other signalling systems including the PhoP/ 
PhoQ system and the RcsB system (21,22). The action of 
the AR regulators is contingent on the growth medium, 
the way in which acid shock is induced, and the growth 
phase, and possibly is affected by strain-specific differ- 
ences. Escherichia coli also shows acid adaptation, chara- 
cterized by enhanced resistance to low pH following 
exposure to mild acidic conditions (23,24). This mechan- 
ism is mediated by the up-regulation of acid shock pro- 
teins, including components of the acid response systems 
described above. Thus its study is used as a way of defining 
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important components of the acid resistance mechanisms 
and their regulation under non-lethal conditions. The 
regulation of the AR2 system has been intensively studied 
in E. coli K-12 and is surprisingly complex, with inputs 
from several different regulatory proteins, including 
global regulators and small RNAs (25). 

Exposure to acid induces a sudden drop at intracellular 
pH (5), which despite the beneficial effect of the ARs, has 
profound effects on the physiology of the cell. For 
example, high intracellular proton concentration may 
induce uncoupling of oxidative phosphorylation resulting 
in alteration of energy metabolism (7). More generally, 
genes involved in energy metabolism, transport and 
amino-acid biosynthesis are known to be modulated at 
both mRNA (26) and protein level (27) suggesting that a 
much broader spectrum of adaptation pathways may be 
modulated in response to acid exposure. Therefore, the 
understanding of acid stress response goes beyond the 
study of the ARs. 

A number of genome-wide expression profiling studies 
representing different acid stress conditions have been re- 
cently published (28-30). Although they have contributed 
to the identification of novel genes transcriptionally 
regulated during acid exposure, they do not yet provide 
us with a comprehensive model of E. coli acid resistance. 

In order to address this challenge, we have taken a 
Systems Biology approach to investigating the gene net- 
works involved in acid resistance. We have integrated 
phenotypic and genome-wide transcriptional data derived 
from wild-type and mutant strains and by using compu- 
tational techniques that allow the inference of regulatory 
networks from observational data (31,32), we have at- 
tempted to identify molecular pathways and regulators 
required for acid adaptation. We show here that failure 
to modulate the expression of metabolic genes, particular- 
ly these encoding enzymes involved in aerobic and anaer- 
obic metabolism, is a dominant feature of strains with a 
strong acid sensitive phenotype. In addition, the applica- 
tion of network inference techniques led us to hypothesize 
that OmpR may be a key regulator of the complex tran- 
scriptional programme involved in acid adaptation. 
Experimental validation of our model, based on the 
analysis of a AompR strain, supported this hypothesis. 

MATERIALS AND METHODS 

Study overview 

The overall aim of this study was to identify pathways 
required for survival in strong acid conditions and to 
discover key regulators of the transcriptional response to 
acid exposure. We first used expression profiling with 
regular sampling over time to characterize the transcrip- 
tional response of a wild-type strain of E. coli BW251 13 to 
mild acid conditions under constant growth conditions 
at an OD 600 = 2. This allowed us to characterize the dy- 
namics of the E. coli response to acid and to identify genes 
differentially regulated over time. By examining the con- 
sequences of deleting 38 of these genes we identified genes 
and pathways that were differentially expressed in mutants 
with different phenotypic strengths. The subsequent 



application of a network inference technique (ARACNE) 
(31,32) allowed the identification of the most likely regu- 
lators of these functions. The hypothesis generated from 
the model was verified by molecular and phenotypic 
analysis of the strain mutated in the putative regulator. 

Bacterial strains 

All the experiments described here were based on E. coli 
K-12 BW25113, which is directly derived from BD792, 
itself a two-step descendent of the E. coli K-12 ancestral 
strain (33,34). All mutant strains analysed in this study 
originated from the Keio collection (35) and were checked 
by PCR to verify the presence of the deletion before being 
used. For this purpose, a combination of locus and 
kanamycin-specific primers were used as described in the 
original publication (36). 

Culture conditions 

Bacterial strains were cultured in Luria-Bertani (LB) 
medium (Sigma Aldrich, USA) supplemented with kana- 
mycin (50 ug/ml). In all our experiments, pH was adjusted 
to neutrality using sterilized 5 M NaOH. In the acid adap- 
tation and acid shock experiments, pH was adjusted using 
sterilized 1 M HC1. Media were buffered with 10% MES. 

In order to maintain optical density within a narrow 
range (1.85-2.15) and to keep growth rate constant 
during the experiment, we used a medium replenishment 
strategy in which medium and cells were removed and 
replaced by an equal volume of pre-warmed medium 
(37°C) at regular intervals (every 5min). Cultures were 
grown for 16 h from a single colony at 37°C in a 
shaking incubator at 200 rpm in 10 ml of LB adjusted to 
at pH 7. An amount of 200 ml of pre-warmed LB medium 
at pH 7 in a 1 1 conical flask was inoculated with overnight 
culture to a starting OD 60 o = 0.1. The culture was grown 
at 37°C and 200 rpm until it reached OD 600 = 2. At this 
point, we initiated the medium replenishment procedure 
by removing 10 ml of culture and adding 10 ml of 
pre-warmed medium at pH 7 every 5 min. Using this pro- 
cedure we kept the culture at OD 600 = 2 and pH 7 for 1 h. 
This ensures that cells reached a steady state before 
addition of acid (Supplementary Figure SI). The pH of 
the culture was then shifted to pH 5.5 by addition of 14 ml 
of 1 M HC1 while maintaining a constant OD for the 
duration of the experiment (1 h). After addition of acid, 
replenishment was performed as before but with buffered 
medium at pH 5.5. 

Phenotypic analysis of mutant strains using flow 
cytometry 

Wild-type and mutant strains were phenotypically 
characterized using flow cytometry analysis (37,38). In 
this application, we used propidium iodide (PI) and 
to-(l,3-diethylthiobarbituric acid) trime thine oxonol 
(BOX) to, respectively, monitor membrane permeability 
and polarization. Briefly, bacterial cells were added to 
2 ml of freshly made, filtered-sterilized (0.2 uM filter) 
PBS supplemented with EDTA, PI and BOX (final 
working concentrations of 4mM, 5 and 10 ug/ml, respect- 
ively). Percentages of healthy and stressed cells were 
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derived from the fluorescence emitted by the cell popula- 
tions after staining, in relation to reference samples of 
alive or dead (ethanol-treated) cells (Supplementary 
Figure S2). The likelihood that the frequency of cells in 
each quadrant of the Flow Cytometry plots was signifi- 
cantly different between the mutants respect and the wild- 
type strains has been computed by using Fisher's exact test 
(39). Contingency tables were constructed using the num- 
ber of healthy cells in the mutant and wild-type strains. 

Plasmid construction 

To generate the complementation plasmid pZCompR the 
complete ompR gene and the native promoter were 
amplified by PCR from the BW25113 chromosome using 
primers ompR-348F (GGTTGCTCGAGCGCCCAGACT 
TGCGGCCCAGG) and ompR + 720R (GGTTGGGATC 
CTCATGCTTTAGAGCCGTCCGG) that introduce 
unique XhoI-BamHI restriction sites. The fragment was 
introduced into the multiple cloning site of the low copy 
number plasmid, pZC320 (40). 

Expression profiling by microarray 

The 10 ml samples of cultures were stabilized by adding 
them to a phenol-ethanol solution (final concentration of 
19% phenol and 1% ethanol). Samples were left on ice for 
20min and then centrifuged (4°C and 2600g) for lOmin to 
recover cell pellets. Stabilized cells were recovered and 
stored at — 80° C. RNA was isolated using the Quiagen 
RNeasy® kit (Quiagen, USA) according to the manufac- 
turer's instructions. Ten micrograms of input RNA was 
labelled using Cy5 labelled dCTP using the CyScribe 
Post-Labelling Kit and purified using CyScribe purifica- 
tion Kit according to the manufacturer's instructions. All 
kits and dyes were from Amersham Biosciences, USA. 
Operon E. coli Ultra GAPS microarray slides (Corning, 
USA) were hybridized overnight with 80 pmol of labelled 
cDNA. The slides were washed in AdvaWash automated 
washing station (Adavlytix, USA) and scanned with the 
ScanArray® GX (PerkinElmer®, USA), using the 
ScanArray® software. 

Quantitative PCR 

To validate gene expression profiles results, E. coli cells 
were grown and stabilized as for the microarray samples. 
RNA extraction was performed as described in the expres- 
sion profiling section. An amount of 40 ng of cDNA were 
analysed with SYBR-Green method, after reverse tran- 
scription of the RNA with the 2000 U Superscript 
II-Invitrogen kit, according to the manufacturer's instruc- 
tions. The primers were designed using Primer 3Plus 
(http://www.bioinformatics.nl/cgi-bin/primer3plus/ 
primer3plus.cgi) and verified for specificity with Primer- 
B LAST (http : / /www . ncbi . nlm . nih . go v/ tools/primer- 
blast/). The sequences for the primers were as follows: 
csr A (left primer CGAGTTGGTGAGACCCTCA, right 
primer AGCCTGGATACGCTGGTAGA); ompR (left 
primer CGTCGCTAATGCAGAACAGA, right primer 
GGTCCACTTCTTCCCCTTTC). Invitrogen primers 
(25nmol) were used at a final concentration of 10 uM. 
SYBR-Green mix from ABGene, with ROX as passive 



reference, was used in a final volume of 10 ul. The 
analysis was performed with a 7900HT Fast Real-Time 
PCR System. Three technical and biological replicates 
were analysed and the gene expression of ompR was con- 
sidered relative to the gene expression of csr A, the expres- 
sion of which does not change under acid conditions 
(41,42). The three biological replicates were then analysed 
with the SDS software (Applied Biosystems) and a Mest 
was performed for significance. 

Data analysis 

The single channel array data were normalized using 
quantile normalization (43) in order to correct for system- 
atic errors. Multivariate exploratory analysis of the gene 
expression data sets was performed using a combination 
of clustering and principal components analysis (PC A). In 
order to identify clusters of genes with similar expression 
profiles, we used SOTA, a clustering methodology based 
on an appropriate distance function (44) with Pearson 
correlation coefficient as a similarity measure. The rela- 
tionships between the transcriptional states of the different 
cell populations were represented using the first two prin- 
cipal components defined by PC A (45). The number of 
components was chosen to represent at least 80% of the 
total sample variance. Bacterial survival and death kinet- 
ics, defined by flow cytometry analysis, were visualized 
using a standard average linkage hierarchical cluster ana- 
lysis (46) with Pearson correlation as a similarity measure. 
Both clustering and PCA were performed using the 
TMEV software application (47,48). 

Differentially expressed genes in the time course experi- 
ment (Figures 1 and 2) were selected using a fold-based 
rule (the absolute value of the log2 of the ratio between 
gene expression intensities between pH 5.5 and pH 7 larger 
than 1.5 at any time point, which corresponds to a 3-fold 
increase or decrease in a linear scale). This identified 2137 
differentially expressed genes. Whenever experimental rep- 
licates were available, differentially expressed genes were 
identified using Significance Analysis for Microarrays 
(SAM) (49) as implemented in the TMEV software appli- 
cation (47,48) by employing a 10% FDR threshold, unless 
otherwise specified in the text. 

Gene lists obtained by cluster or differential expression 
analysis were assessed for over-representation of KEGG 
(50) and Gene Ontology (51) functional terms using the 
open-source software DAVID (52,53) (http://david.abcc 
.ncifcrf.gov/). The GO analysis was performed with the 
lowest level of GO terms. In all cases, we used an 
FDR < 1 % threshold to define a statistically significant 
enrichment. 

Gene networks representing the neighbourhood of the 
two-component system regulators were inferred using the 
software application ARACNE (31,32) using the expres- 
sion matrix representing the transcriptional state of 27 
mutant strains and three replicates of the wild-type 
strain profiled at pH 7 and 5.5 (for a total of 60 arrays). 
ARACNE is an algorithm used for the identification of 
transcriptional interactions between gene products, by 
identifying statistical interactions based on the mutual in- 
formation (32). Significant interactions were defined by a 
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Figure 1. Principle component analysis of the transcriptional response of BW25113 to acid adaptation. This analysis shows the change in the 
transcriptional response of BW25113 in the first hour of acid adaptation. The x and y axes represent, respectively, the first and second components 
in the PCA. 



P-value threshold of P < 10 -7 (corresponding to one over 
2000 false positive connections). In order to eliminate 
non-direct interactions, we used the inequality principle 
as implemented in ARACNE with a DPI of 0.1 (32). 
The resulting networks were visualized using the 
Cytoscape software application (54). 

To identify potential binding sites for OmpR, we have 
used the online tool Virtual Footprint (http://prodoric 
.tu-bs.de/vfp/) using the Prodoric library and performing 
the Regulon and the Promoter analysis (55). The PWM 
(Position Weight Matrix) considered were OmpC-box and 
OmpF-box of E. coli strain K-12. 

Microarray data 

Gene expression profiles data are accessible on the Gene 
Expression Omnibus (GEO) website (http://www.ncbi 
.nlm.nih.gov/geo/), with accession number GSE13361. 

RESULTS 

Escherichia coli adaptation to acid involves a rapid 
but transient transcriptional response 

The approach we took in this study was to first character- 
ize the dynamics of the transcriptional response of E. coli 
to acid adaptation. In order to do this, we monitored 
cultures that were kept in a balanced state of growth by 



removal of culture and replenishment with pre-warmed 
medium at regular intervals, a procedure that maintains 
cells in a constant transcriptional state (Supplementary 
Figure SI). The effect of acidification on gene expression 
was analysed using microarrays with RNA samples from 
six time points: prior to treatment, and 30 s, 5, 15, 25 and 
60min after exposure to pH 5.5. A high-level representa- 
tion of the changes in the cells' molecular state, performed 
using PCA, revealed that the process is defined by a rapid 
(30 s to 15 min) but transient response leading, after 1 h, to 
cells whose transcriptional state is similar to cells grown at 
pH 7 (Figure 1). 

A more detailed analysis of the gene expression dynam- 
ics both by visual inspection of the data and using clus- 
tering revealed that the transcriptional events during this 
period of adaptation could be summarized by four 
main clusters of gene expression profiles (Figure 2). The 
first two clusters represented the earliest response to 
stimulation, including down- and up-regulated genes 
(Figure 2, clusters 1 and 2, respectively). In these two 
clusters, the largest changes in absolute value were 
detected 5 min after acid exposure. The second set of 
two clusters represented genes that are more gradually 
modulated down or up in response to pH 5.5 under 
these conditions (Figure 2, clusters 3 and 4, respectively) 
and where the largest change in absolute value was seen 
25 min post-treatment. 
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Figure 2. Cluster analysis of gene expression profiles in response to acid adaptation. This Figure shows the result of cluster analysis performed using 
SOTA (41) on the most differentially expressed genes. SOTA identified four main clusters of co-ordinately regulated genes. A heatmap representing 
the expression profile of all genes in each cluster is represented on the left side of the figure (green and red correspond respectively to expression levels 
below and above the mean of gene expression across the time points). Clusters 1 and 3 show genes that are transiently down-regulated whereas 
Clusters 2 and 4 show genes that are transiently up-regulated. The table to the right of the heat maps shows GO and KEGG functional terms 
significantly enriched in each cluster (Count, number of genes for each pathway; Benjamini, P-value correction; ASRs, Acid Stress Response systems 
genes; AFI, Acid Fitness Islands genes). 



We looked to see whether genes in each cluster were 
significantly enriched for any specific GO and KEGG 
functional terms (Figure 2). Genes represented in the sig- 
nificant terms were then mapped on the BW251 13 genome 
and, as an additional level of quality control, the tran- 
scriptional response of genes in operons was analysed to 
check that all genes in a given operon showed that same 
transcriptional response (Supplementary Figures S3A-D 
and Figure S4: 17). 

Consistent with our current understanding of acid 
adaptation, we observed the up-regulation of many of 
the genes associated with the glutamate, arginine and 
lysine-dependent acid response systems, and members of 
the AFI such as gadW, hdeA, hdeD, hdeB, yhiD, sip and 
yhiF. The exceptions found in this preliminary analysis 
were the genes coding for the isozyme GadA, and the 
GadX and the GadE regulators, since during early sta- 
tionary phase they are supposed to be active, following 
the expression of RpoS (56). This apparent inconsistency 
may reflect a more complex regulatory system driving 
the transcriptional behaviour of the GAD genes. Genes 
coding for the FoFl ATPase complex (thought to be a 
component of the amino acid-independent acid stress 
response system AR1) were also down-regulated. In 
addition, the transcriptional response to acid exposure 
involved the concomitant up-regulation of genes involved 
in anaerobic respiration (GO:0009061) and down- 
regulation of genes involved in aerobic respiration 
(GO: 0009060). Genes involved in the regulation of the 
cell wall (GO:0005618) and translation (GO:0006412) 
were down-regulated whereas genes involved in membrane 
transport (ecd02010) were up-regulated in response to acid 
exposure (Table 1 and Figure 2). We also found that 
several genes coding for membrane proteins involved in 
osmoregulation were up-regulated (Supplementary 
Figures S6 and S10). More specifically, a number of 
osmoprotectant transporters involved in the response to 
hyperosmotic shock were transiently up-regulated. More 
specifically, a number of osmoprotectant transporters 
involved in the response to hyperosmotic shock were tran- 
siently up-regulated together with the genes mscL and 
mscS coding for the mechano-sensors involved in the 
response to hypo-osmotic shock. 

As this initial analysis of the time course data was car- 
ried out using a single experiment only at each time point, 
we also compared three replicates of wild-type control 
cells with cells 15min after acid exposure. Statistical 
analysis of these data identified 1871 differentially ex- 
pressed genes. The result was largely overlapping with 
the results of the time course analysis (83% of the genes 
identified in the time course experiment were also ident- 
ified by the single time point analysis, Supplementary 
Figures S5 and S6), showing our time course data to be 
robust. 

Genes transcriptionally regulated by mild acid 
are generally required for surviving acid shock 

Having performed an initial characterization of the E. coli 
transcriptional response during acid adaptation, we 
wished to see what fraction of the differentially expressed 
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Table 1. Gene expression in response to acid adaptation 



Function 


Complex 


Genes 


Regulation 


P- value 


• 

Oxidative 


NADH dehydrogenase 


nuoA 


I 


1.18E-1 1 


phosphorylation 




nuoB 


I 








nuoC 


I 








nuoE 


I 








nuoF 


I 








nuoK 


I 






tori A 1 rase 


atpB 


I 








atpE 


I 








atpH 


I 








atpA 


I 




Glycolysis 


Galactose- 1 -epimerase 


galE 


I 


i cor? A A 

3.58E-04 






galT 


I 








galK 


I 








galM 


I 






Phosphoglycerate kinase 


epd 


I 








Pgk 


I 




Translation 


Glycil-tRNA synthetase 


glyS 


I 


1 f\1U AC 

3.U/E-ID 






glyQ 


I 






Lysyl-tRNA synthetase 


lysS 


I 








dsbC 


I 




Anaerobic 


Formate dehydrogenase 


fdnG 


t 


L.20h-0o 


respiration 




fdnH 


t 








fdnl 


t 






Trimethylamine N-oxide 


torC 


t 






reductase 


tor A 


t 








torD 


t 




Transport 


Iron dicitrate ABC 


fecA 


t 


1.20E-06 




transporter 


fecB 


t 








fecC 


t 








fecE 


t 






Sulphonate/nitrate/ 


ssuA 


t 






taurine ABC 


ssuB 


t 






transporter 


ssuC 


t 








ssuD 


t 






Probable proton-driven 


yjeP 


t 






drug efflux system 


hsrA 


t 





The table list genes belonging to some of the pathways that are down 
(oxidative phosphorylation, glycolysis and translation) and 
up-regulated (anaerobic respiration and transport). The P-values 
obtained from DAVID functional annotation analysis are shown in 
the right side of the table alongside the direction of change. 



genes was required for survival during acid shock. In order 
to address this question, we selected 38 genes representa- 
tive of the main functions modulated during acid adapta- 
tion (Supplementary Table SI) and tested the ability of 
strains where each of these genes was deleted to survive 
in strong acid conditions (pH 2.5). The mutants were 
assayed at different time points (5, 30, 60, 90, 120, 150 
and 360 min) after acid exposure, and alive and dead 
cells were quantified by flow cytometry analysis. All 
tested mutants showed a detectable increase in sensitivity 
to low pH in the assay, though some effects were small 
(Figure 3, Supplementary Tables S2 and S3). 

In order to define the severity of the phenotype, we 
chose a threshold, which allowed us to cluster the 
mutants in three different groups as having strong, weak 
or intermediate phenotypes. We considered 30 min of ex- 
posure at pH 2.5 the parameter for the definition of the 
phenotype: strong phenotype mutants show < 3 5-40% of 
healthy cells after 30 min at pH 2.5, intermediate pheno- 
type show between 60% and 40% of healthy cells, while 
weak phenotype show >60% of healthy cells (Figure 3). 
The wild-type strain, included for comparison, shows the 
weakest phenotype, as expected, as this strain is quite acid 
resistant under the growth conditions of this assay. 
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Figure 3. Clustering of mutant strains according to their phenotypic response to acid shock. The figure represented in the format of a heatmap the 
percentage of healthy cells following exposure to LB at pH 2.5, without prior adaptation for the mutants we tested. Survival was measured at 
different time points (one sample every 30min for 3h). The colour code indicating the percentage of survival determined by flow cytometry analysis 
at each time point is below the heatmap. Hierarchical clustering of the survival time course data revealed three clusters of strains on the basis of their 
survival profiles. These were labelled as weak, intermediate and strong phenotypes. 



Genes that when mutated gave a strong phenotype 
include transcriptional regulators (phoP, rpoS and ihfB), 
a glycine decarboxylase (gcvP), the 2-oxoglutarate de- 
hydrogenase (sucB) and the two-component system sensor 
phoQ. Strains with an intermediate phenotype included 
mutants of two acid stress chaperones (hdeA and hdeB) 
and nitrate reduction (nar J and narX). Strains with a weak 
phenotype included mutants of several transcriptional 
regulators (crp, rpoD, hns,phoB, ihfA,fliA,fnr) and, inter- 
estingly, included three of the GAD system-specific regu- 
lators (gadE, gadW and gadX). 

Defective expression of cell wall and energy metabolism 
genes correlates with increased sensitivity to acid 

The experiments described above showed that the inacti- 
vation of genes that are transcriptionally regulated during 
adaptation often results in a significant loss of survival 
following acid shock. We reasoned that analysis of the 



transcriptional response to acid exposure in the different 
phenotypic groups would enable us to formulate hypo- 
theses concerning the molecular pathways that are more 
important for survival of acid stress. We, therefore, sub- 
jected the same 38 mutants used above to a microarray 
analysis after acid adaptation and asked whether any com- 
ponent of the transcriptional response was predictive of 
loss of acid resistance. For practical purposes, this experi- 
ment was performed on a single time point (15min after 
exposure at pH 5.5), which is where the highest change in 
gene expression in the early response clusters of the wild- 
type strain occurred (Figure 2). A PC A analysis confirmed 
that the acid adaptation transcriptional programme was 
defective in the mutant strains (Supplementary Figure S8). 

In order to identify molecular pathways linked to the 
severity of the phenotype, we compared the gene expres- 
sion profiles (expressed as a log-ratio between expression 
values at pH 5.5 and 7) between the wild-type (represented 
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by four replicated experiments) and mutated strains that 
showed either an intermediate or a strong phenotype. We 
found 221 genes to be differentially expressed between the 
three experimental groups (FDR< 10%). Among the list 
of genes up-regulated in the wild- type strain, we found 
some of them involved in anaerobic respiration (hyfC, 
nrfF, men A), sugar transport (xylE, ulaA, ycjP, sgcC, 
malF, agaD, agaV), lipopolysaccharide biosynthesis (ylbH, 
rfaF, rfe, rfaJ, rhsA) and purine metabolism (yqeA, nudF, 
allB), whereas the list of down-regulated genes was en- 
riched in aerobic respiration genes (including oxidative 
phosphorylation and TCA cycle) (ubiA, nuoF, acnB, 
sdhC, acnA, sdhD, ppc, nuoN, cyoA, cyoB, cyoD), DNA 
repair (uvrC, xthA, nei, exoX) and cell cycle (mrdB, ftsY, 
mukE, ftsA) (Figure 4). 

We observed that the large majority of the genes pre- 
dictive of phenotypic outcome (89%) were also 

WT INT STR 




• Anaerobic respiration 

- • Carbohydrate transport 

• Cell envelope biosynthesis 

• Purine metabolism 



n 




• Aerobic respiration 

• Oxidative phosphorylation 

• TCA cycle 

• Cell cycle 



log 2 Ratio 

Figure 4. Differential gene expression in response to acid in the three 
phenotypic groups. The figure represents in a heatmap format the ratio 
of gene expression in acid-exposed cells respect to control cultures of 
the 221 differentially expressed genes between the wild-type and the 
mutants belonging to the intermediate (INT) and strong (STR) pheno- 
type groups. The annotation on the right refers to functional groups of 
genes that were either down regulated (green) or up-regulated (red) in 
the wild-type strains at pH 5.5 relative to their expression at pH 7. 
These are shown in the first column of the heat map. In the second and 
third columns, the means of the log ratios of expression at pH 7 and 
5.5 for these same genes are shown for the intermediate and strong 
groups. 



differentially expressed in response to acid exposure 
(Figure 4). Remarkably, the ratios of the expression be- 
tween cells grown at pH 5.5 and 7 were generally reversed 
in mutant strains that had an intermediate phenotype, 
while mutant strains with a strong phenotype were 
characterized by a general inability to regulate most of 
these genes (Figure 4). 

Network inference analysis reveals new potential 
regulators of the global transcriptional response 
to acid stress 

The results above identified molecular pathways whose 
levels of expression are linked to survival during acid 
shock. We therefore investigated which transcription 
factors might regulate the activity of these pathways. In 
order to focus on regulators whose activity may be directly 
influenced by environmental sensing we specifically looked 
for direct or indirect targets of the regulatory subunits of 
two-component systems. To achieve this, we used an 
unbiased approach that relies on a well-validated reverse 
engineering methodology (ARACNE). ARACNE uses a 
probabilistic measure of dependency (mutual information) 
to identify potential targets of a given transcription factor. 
It relies entirely on the analysis of a compendium of ex- 
pression profiles representing different perturbations of a 
cell's homoeostatic state. Genes are inferred to be tran- 
scriptionally coupled when there is a statistically signifi- 
cant correlation between their levels of expression, across 
all samples in the compendium data set and after potential 
indirect connections are removed using the inequality 
principal criteria (32). Typically, large compendia of 
microarray data are used for the inference procedure. In 
this analysis, we chose to use the collection of 60 experi- 
ments representing expression profiles of genes in mutants 
and wild-type strains at pH 7 and 5.5, described in the 
previous paragraph. This had the advantage that it repre- 
sented a collection of mutant strains relevant for acid 
adaptation while being sufficiently large to allow reliable 
inference (57). Since our objective was to identify regula- 
tors of the previously identified pathways, we ranked all 
the E. coli two-component system regulators on the basis 
of the number of connections they had with genes belong- 
ing to those pathways that we had shown to be predictive 
of phenotypic response and highly regulated during acid 
response. The ranking is shown in Table 2. The highest 
ranking gene by this analysis was OmpR, the regulatory 
subunit of the osmoregulator two-component system 
EnvZ/OmpR. OmpR was connected to genes involved 
in aerobic energy metabolism (pyruvate metabolism and 
glycolysis), signal transduction and transport, as well as 
some of the components of the GAD system (Figure 5 
and Supplementary Table S5). NarP, the nitrate/nitrite 
response regulator, was connected to nucleotide and 
amino acid metabolism; the envelope stress response regu- 
lator BaeR was connected to genes involved in transport 
and metabolism. Most of the other two-component system 
regulatory subunits were poorly connected to genes in the 
pathways we had previously identified. On the basis of 
these observations, we proposed the novel hypothesis 
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Table 2. Potential two-component systems regulators targets as defined by the ARACNE analysis 



Gene 



Two-component system 



Targets 
based 
on MI 



Function enriched 



P-value 



ompR 



Omsoregulatory two-component system 
OmpR/EnvZ 



33 



dcuR 

baeR 

narL 
phoP 
kdpE 

rstA 

uvrY 

narP 

rcsB 



fimZ 
torR 



rssB 

yfjR 

uhpA 



Two-component system DcuR/DcuS, 
regulating anaerobic fumarate 
respiratory system 

Two-component regulatory system 
BaeS/BaeR 



NarX/NarL Two-component system, 

nitrate/nitrite dependent 
PhoP/PhoQ Two-component system, 

magnesium dependent 
Two-component regulatory system 

KdpD/KdpE involved in the 

regulation of the kdp operon 
Two-component regulatory system 

RstA/RstB 

Two-component regulatory system 
UvrY/BarA 

NarP/NarQ Two-component system, 
nitrate/nitrite dependent 



RcsB/RcsC regulates the expression of 

genes involved in colanic acid capsule 

synthesis 
involved in fimbrial expression 
TorR/TorS two-component system, 

responding to changes in the 

concentration of the TMAO 
facilitates and regulates degradation 

of sigma S 
Two-component system YfjR/YfhK 
UhpA/UhpB two-component system 

responding to external concentrations 

of glucose-6-phosphate 



27 

43 

71 
48 
31 

27 

23 

50 

36 



34 
20 



18 



13 



GO_BP intracellular pH elevation (3 genes) 
GO_CC cytoplasm (7 genes) 
KEGG ABC transporters (3 genes) 
KEGG gluthatione metabolism (2 genes) 
KEGG glycolysis/gluconeogenesis(4 genes) 
KEGG pyrimidine metabolism (5 genes) 
KEGG pyruvate metabolism (2 genes) 
KEGG two-component system (3 genes) 
GO_BP ciliary or flagellar motility (3 genes) 
KEGG two-component system (7 genes) 
KEGG methane metabolism (3 genes) 
KEGG pyruvate metabolism (3 genes) 
GO_BP transport (10 genes) 
KEGG ABC transporters (11 genes) 
KEGG galactose metabolism (5 genes) 
KEGG purine metabolism (3 genes) 
GO_MF metal ion binding (4 genes) 
KEGG glycolysis/gluconeogenesis (3 genes) 
GO_BP oxidation reduction (3 genes) 
KEGG oxydative phosphorylation (3 genes) 
GO_BP nickel ion transport (2 genes) 
KEGG two-component system (6 genes) 

GO_BP translation (4 genes) 

KEGG alanine, aspartate and glutamate metabolism (3 genes) 
KEGG Ribosome (3 genes) 

GO_BP two-component signal transduction system (5 genes) 

GO_BP transcription (7 genes) 

KEGG two-component system (3 genes) 

GO_BP pyrimidine nucleotide biosynthetic process (6 genes) 

GO_BP iron-sulphur cluster assembly (6 genes) 

KEGG pyrimidine metabolism (7 genes) 

KEGG alanine, aspartate and glutamate metabolism (5 genes) 
GO_BP transcription (9 genes) 



GO_BP transcription (7 genes) 
GO_BP metabolic process (5 genes) 

GO_MF protein binding (8 genes) 

GO_BP transcription (5 genes) 
GO_CC integral to membrane (8 genes) 



4.50E-05 
6.80E-03 
4.20E-02 
3.90E-02 
3.00E-03 
3.60E-03 
2.40E-03 
3.10E-03 
1.10E-03 
3.10E-05 
9.28E-05 
3.46E-03 
1.25E-06 
3.21E-06 
1.39E-05 
2.06E-03 
2.88E-03 
4.63E-04 
1.63E-03 
1.63E-03 
5.78E-04 
4.48E-04 

5.36E-03 
9.81E-04 
4.27E-05 
2.48E-05 
1.45E-04 
4.27E-05 
3.36E-10 
1.21E-07 
1.72E-06 
1.57E-05 
1.08E-04 



7.24E-04 
1.09E-02 



7.77E-03 

9.98E-03 
4.33E-03 



Regulators were ranked on the basis of how many connections were found to functions that were modulated during the acid adaptation time course 
(in red). P-values of the DAVID functional annotation analysis are shown. 



that OmpR may be a key regulator of acid response in 
E. coli BW25113 under the conditions of our study. 

Experimental confirmation that OmpR is a key regulator 
in the response to low pH 

If our hypothesis is correct, we expect a mutant strain 
lacking ompR to display the following properties. First, 
the genes differentially regulated between AompR and 
the wild-type strain should significantly overlap with 
genes differentially regulated during acid adaptation in 
the wild-type. Second, this overlap should be consistent 



with the regulation of ompR itself during acid adaptation 
in the wild-type strain. For example, if ompR expression 
decreases on acid shock, then those genes that are 
normally repressed by OmpR should increase after acid 
shock, and show significant overlap with genes that are 
up-regulated if ompR is deleted. Third, the AompR mutant 
strain should be less able to initiate the normal transcrip- 
tional response to acid adaptation and, hence, should 
show a significant decrease in survival of acid shock 
relative to the wild-type strain. In order to test these pre- 
dictions, we analysed the expression profile of a AompR 
mutant strain at pH 7 and 5.5, and performed a 
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Figure 5. Network of regulatory interactions in the neighbourhood of OmpR, inferred using ARACNE. Nodes represent genes and edges represent 
inferred gene-to-gene connections. Genes are colour coded on the basis of their function. 



phenotypic characterization of this strain by flow 
cytometry after a direct challenge at pH 2.5. The results 
of this analysis closely matched our predictions. Genes 
differentially regulated between wild-type and the 
AompR mutant strain showed significant overlaps with 
genes differentially expressed in the wild-type strain 
during acid adaptation (Figure 6C and D). More specific- 
ally, 120 (23%) of the up-regulated genes in the wild-type 
during acid adaptation were also expressed more highly in 
the AompR mutant than in the wild-type at pH 7, while 
280 genes (71%) genes down-regulated in the wild-type 
during acid adaptation were expressed more weakly in 
the AompR mutant than the wild-type at pH 7. The dir- 
ection of the overlap was consistent with the observed 
down-regulation of ompR in response to acid exposure 
(Figure 6 A and B). qPCR showed the down-regulation 
of ompR at pH 5.5 with values of slope and r 2 for the 
efficiency and the accuracy of the measurements: csrA 
slope = —3.1 cycles/logdecade, r 2 = 0.99; ompR slope = 
—3.01 cycles/logdecade, r 2 = 0.99. Even more strikingly, 
there was no transcriptional response detectable to acid 
exposure in the AompR mutant either using a univariate 



statistical analysis approach (no differentially expressed 
genes were identified up to a FDR < 50% threshold) or 
using a multivariate exploratory analysis approach (Figure 
6D). Flow cytometry analysis showed that the AompR 
mutant strain showed one of the strongest phenotypes 
detected in our study (Figure 7, Panel A). Mutant cells 
showed increased resistance to acid shock if pre-adapted 
to mild acid conditions, suggesting a small residual ability 
to mount an effective response, although they were still 
significantly less resistant than the unadapted wild-type 
strain. Complementation of the AompR mutant strain 
with a copy of the wild-type ompR gene expressed under 
its own promoter in the low-copy number plasmid pZC320 
led to restoration of normal levels of acid resistance, 
whereas the vector alone had no effect [Figure 7, Panel 
B shows analysis by flow cytometry; all results were also 
confirmed with plate counts (data not shown)]. 

Comparison between effects of acid exposure in different 
mutant strains confirms the importance of OmpR 

To compare the effects of different mutations on the acid 
resistant phenotype with that of the AompR mutation, we 
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Figure 6. Comparison of transcriptional responses of AompR and wild-type to acid stress of pH 5.5. Panels (A and B) show the expression levels of 
ompR at pH 7 and after 15min exposure at pH 5.5, based on either microarray data (A) or qPCR data (B). Bars show standard deviations of three 
biological replicates. ompR is significantly down-regulated at pH 5.5 in both data sets (P-value 0.02 and 0.01, respectively, for microarrays and qPCR 
data, obtained with /-test). (C and D) show the extent of overlap between genes which are differentially up-regulated (C) or down-regulated (D) in 
the wild-type at pH 5.5, relative to expression at pH 7, and genes which are over-expressed (C) or under-expressed (D) in AompR at pH 7, relative to 
the wild-type at pH 7. (E) PCA plot of transcriptome changes in wild-type and the AompR mutant, analysed at pH 7 and 5.5. Blue dots, wild-type 
pH 7, pink dots: wild-type pH 5.5, black dots AompR at pH 7, red and orange dots AompR at pH 5.5 after 30 s and 15min of exposure, respectively. 



did flow cytometry experiments on BW25113 strains 
deleted for gadE, gadC, adiC and rpoS. GadE is a 
central regulator of the GAD system, GadC is the 
GABA/glutamate antiporter needed for function of the 
AR2 system, AdiC is the arginine/agmatine antiporter 
needed for function of the AR3 system, and RpoS is a 



sigma factor of RNA polymerase involved in the acid re- 
sistant phenotype. Under our conditions, we confirmed 
that the /AompR strain showed the strongest phenotype 
of all these mutants when cells were not adapted with a 
prior shock at pH 5.5 (Figure 7, Panel C). After lOmin of 
exposure, <30% of the cell in /AompR were still alive, 
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Figure 7. (A) Acid sensitive phenotype of the AompR strain. Survival of wild type and AompR strains in response to exposure to pH 2.5, either 
directly or after prior adaptation to pH 5.5, was measured by flow cytometry. Error bars show standard deviations of four independent biological 
replicates. The y axis shows % survival relative to viable cell numbers at t 0 . Wild type without adaptation, pale blue; wild-type with adaptation, dark 
blue; AompR without adaptation, red; AompR with adaptation, brown. (B) Complementation of the AompR strain restores acid resistance. The 
AompR mutant was complemented with the plasmid pZCompR, or as a control with the empty vector pZC320, and acid resistance measured at pH 
2.5 without prior adaptation Error bars show standard deviation values of four biological replicates for WT and AompR and three biological 
replicates for the complemented strain. Wild type, pale blue; AompR without vector, red; AompR with empty vector, green, AompR with pZCompR, 
mauve. (C, D) Comparison of the AompR acid resistance phenotype with those of other mutations in key acid response genes. Survival curves of the 
following mutants after exposure to pH 2.5, without (A) or with (B) prior adaptation at pH 5.5 for one hour: AompR (red), AgadE (yellow), AgadC 
(purple), AadiC (green) and ArpoS (dark blue) compared to the wild type (blue). Error bars show standard deviation values of four independent 
biological replicates for wild-type and AompR and three biological replicates for the other mutant strains. 



while in the other strains, >70% of the cells were healthy 
at this stage. Strains carrying mutations in either of the 
two amino acid antiporters showed a weak phenotype. 
Different effects were seen when cells were induced 
at pH 5.5, before exposure to acid shock at pH 2.5 
(Figure 7, Panel D). All mutants in this case (apart from 
the adiC mutant) showed large reductions in viability 
compared to the wild-type strains; however, the AompR 
strain still had the most severe phenotype. 

DISCUSSION 

Our work has shown that the ability of E. coli BW251 13 to 
survive in strong acid conditions is dependent on an 
OmpR-dependent transcriptional programme that 
involves the regulation of energy metabolism. Mutant 
strains affecting the classic acid resistance genes generally 
have a weaker phenotype than a /AompR strain, suggesting 



that this novel response programme is the key to acid 
adaptation. 

A rapid shift between the expression of genes involved in 
aerobic and anaerobic energy metabolism is a key 
landmark of acid adaptation 

Several lines of evidence that emerge from this work sup- 
port the view that acid exposure induces a shift between 
aerobic and anaerobic metabolism and that this is a strong 
requirement for survival during acid shock. First of all, 
E. coli BW25113 cells under the conditions of our experi- 
ment express high levels of mRNA for enzymes involved 
in aerobic metabolism and lower levels of genes involved 
in anaerobic energy metabolism. A shift to low pH cul- 
ture conditions induces a reversal of this pattern 
(Supplementary Figure S4: 1-7). The functional analysis 
revealed a trend in the modulation of the genes belonging 
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to the anaerobic respiration pathway (according to GO 
terms), which involves the down-regulation of the 
fumarate-dependant respiration and the up-regulation of 
nitrate and formate respiratory enzymes (Table 3). 
Enterobacteria are predominantly fermentative anaerobes 
that generate acids. The dominant metabolism of E. coli in 
an anaerobic environment is based on the mixed acid 
fermentation, and on the adaptation to readily decrease 
external pH, for example, by inducing the synthesis of 
formate hydrogen lyase pathway. The metabolic switches 
imply that acid triggers a shift from aerobic respiration 
driven by the tricarboxylic acid cycle to adaptation to 
oxygen-limited growth. This suggests that E. coli has 
evolved the ability to prepare for anaerobic lifestyle by a 
pre-emptive induction of some anaerobic pathways prior 
to oxygen starvation. The comparison between the 
wild-type and mutant strains (with intermediate and 
severe phenotypes) revealed that failure to invert the ex- 
pression of aerobic and anaerobic metabolism genes cor- 
relates with a strongly acid sensitive phenotype. 
Moreover, many of the strains which show a strong 
phenotype are mutated in genes involved in energy metab- 
olism (Figure 4 and Supplementary Figure S9: 1-2). This 
indicates that under the conditions described here, effect- 
ive proton scavenging alone may not be sufficient for 
survival unless the expression of bioenergetics genes is 
also modulated during adaptation to low pH conditions. 
We have seen that the genes transcriptionally controlled 
by the pH change are essential for the entire process pre- 
viously described. This implies that the gene products 
are part of an integrated adaptive process in which tran- 
scription regulation depends on the completion of the 
process itself. There are currently few insights into the 
mechanism underlying this response. Suggestions can be 
made if an fnr mutant may be defective in adaptation to 
low pH. Although FNR is an iron-sulphur protein that is 
inactivated by oxygen (58,59), in cultures of moderately 
high density such as those used in the current experiments, 
FNR is partially active despite vigorous aeration. It would 
also be interesting to see whether FNR is critical for adap- 
tation to acid of cultures at a much lower bacterial density. 

OmpR: an important acid response regulator 

Our study provides strong evidence for a role of OmpR as 
a regulator of the transcriptional response to acid adapta- 
tion in BW25113. Although OmpR has not been directly 
implicated hitherto in acid resistance in E. coli, it has been 
shown than an ompR UPEC mutant shows reduced 
survival in the mouse urinary tract, and that the growth 
defect seen in this mutant in high salt is enhanced at low 
pH (5.5); lethal pH was not tested in this experiment (60). 
OmpR has been shown to regulate the stationary phase 
acid inducible response in Salmonella Typhimurium (61) 
potentially by counteracting H-NS-mediated repression 
(61,62), and an OmpR-like regulator (HP0166) has also 
been implicated in the acid response of H. pylori (63) sug- 
gesting a broader role for regulators of this type in acid 
stress responses. 

Gene expression in the ompR mutant partially mimics 
the response to acid of the parent strain (Figure 6) but 



mutant cells are unable to mount any response to acid, 
suggesting that modulation of OmpR is required for an 
effective response to acid. This result, coupled with the 
acid sensitivity of the AompR mutant, are consistent 
with a model where OmpR is required for the expression 
of cellular components, or for the establishment of a par- 
ticular cellular state, which is needed for the cells to be 
able to respond to acid stress in a way that enhances their 
survival. In the absence of OmpR, this state no longer 
exists and so the AompR mutant fails to respond to acid- 
ification and shows enhanced acid sensitivity. Cells lacking 
this important regulator are even less able to survive in 
extreme acid conditions than those carrying mutations in 
other important genes implicated in acid resistance. The 
key role of OmpR in regulating adaptation to low pH was 
unexpected and has not featured in an extensive literature 
on the EnvZ-OmpR two-component system. Key ques- 
tions for future research to answer are whether EnvZ 
mediates the response and, if so, to what chemical signal 
EnvZ responds. In a previous study on Shigella, FNR, the 
dual transcriptional regulator of the switch between 
aerobic and anaerobic metabolism, was found to 
regulate the length of type III secretion system needles, 
required for secretion of the invasion plasmid antigen. 
Exposure to oxygen at the surface of the gastrointestinal 
mucosa inactivates FNR, reversing the block on invasion 
antigen secretion and hence priming the bacteria for the 
attack (58,59). Perhaps the OmpR-mediated response to 
acid fulfils a similar role. Moreover, a deeper analysis of 
FNR targets revealed a down-regulation of the functions 
modulated by FNR, which we found important for our 
work (Supplementary Table S6). 

Whether OmpR exerts its function directly or by 
modulating the activity of other regulators remains an 
important question. We found that some of the genes 
downstream from OmpR (as defined by correlation and 
KO analysis) have OmpR-binding sites in their promoter 
regions (Supplementary Figure S7), suggesting that at 
least in some cases OmpR may directly activate genes 
involved in acid response. 

The relevance of the osmotic stress response to acid 
adaptation 

The novel role we have found for OmpR suggests a link 
between osmotic and acid stress responses. A potential 
connection between the activation of anaerobic metabol- 
ism (an important feature of acid adaptation in this strain) 
and osmotic stress has been previously described (64), 
where it was linked to changes in DNA topology. Many 
genes involved in the response to osmotic shock are also 
modulated in response to acid in our experiments 
(Supplementary Figures S6 and S10). Most of the channels 
involved in the transport and production of osmopro- 
tectants were down-regulated at pH 5.5. However, the 
H + /proline symporter (65), the choline transporter (66) 
and two K + channels (67,68) (Supplementary Figure S10) 
were up-regulated after acid exposure, a response that 
makes adaptive sense in light of the finding that K + and 
proline have a beneficial effect on pH homoeostasis in 
E. coli (69). It may also be the case that the import of 
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Table 3. Regulation of the genes involved in the anaerobic respiration pathway, according to GO terms 
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The list includes genes down-regulated in response to acid exposure (left) and up-regulated (right). 



amino acids into the cytoplasm, as part of the acid adap- 
tation response, could drive the cells to a weak condition 
of hypo-osmolarity. Consistent with this hypothesis, the 
mechanosensor channels, MscS (70) and MscL (71), with 
the acquaporine AqpZ (72) and the K + mechanosensor 
(73) were up-regulated in the first 15min after acid expos- 
ure. Substantial changes are known to take place in the 
balance of ions across the membrane on acid shock in 



E. coli, leading to reversal of membrane polarity, and in 
this context it is not surprising that a regulator that itself 
responds indirectly to changed osmolarity is crucial in 
regulating this response. 

A major two-component system involved in regulating 
the acid stress response, the EvgAS system, did not emerge 
from this study. This is unsurprising as it is known that the 
expression of the evgAS operon is little affected by acid 
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and so it would not have been detected as being important 
using the methods described here. Furthermore, the evgAS 
system is mainly thought to be important in regulating the 
gad system, but this system has a limited role in the con- 
ditions described here, as can be seen by the fact that a 
gadE deletion only confers a mild phenotype. 

Different E. coli strains show a wide range of resistance 
to acid. It is not yet clear whether OmpR expression can 
explain some of the variation in this naturally occurring 
resistance. It may be that the relative importance of 
OmpR may depend on the strain being used, and that 
BW25133 is particularly dependent on OmpR. We are in 
the process of testing this hypothesis. 

Overall, three aspects of the physiological response of E. 
coli K-12 to mild acid stress require further analysis and 
explanation. First, the response regulator OmpR, which 
normally is thought of as responding to osmotic stress, is 
essential for adaptation to low pH and for priming 
survival of more severe acid stress. Second, many genes 
previously identified as components of OmpR regulon as 
well as additional previously unrecognized members, are 
essential for adaptation. Third, mild acid stress induces 
various metabolic switches, for example, from glycolysis 
to gluconeogenesis and fatty acid synthesis, and from 
energy generation powered by the TCA cycle to expression 
of FNR-regulated genes associated with anaerobic respir- 
ation (Supplementary Table S6, enriched functions of 
FNR target genes in response to acid exposure in our 
time course data). 



CONCLUSIONS 

We have taken a broadly applicable, systems-level 
approach to defining the response to acid of E. coli 
BW25133. This approach has revealed several entirely 
novel features of this well-studied response, including the 
importance of a shift in the metabolic state of the cell, and 
a central role for the OmpR regulator, a regulator that 
was already know to play an important role in modulating 
E. coli gene expression under standard conditions (74). 
The next steps in this study will be to see how broadly 
applicable these findings are across different strains, 
including pathogens and different growth conditions, to 
dissect the molecular and physiological mechanisms that 
underlie the changes we have detected, and to understand 
in more detail their adaptive value for the cell. 
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